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Abstract 


Large-scale winter teleconnection of the East Atlantic / West Russia (EA/WR) over the 
Atlantic and surrounding regions is examined in order to quantify its impacts on temperature and 
precipitation and identify the physical mechanisms responsible for its existence. A rotated 
empirical orthogonal function (REOF) analysis of the upper-tropospheric monthly height field 
captures successfully the EA/WR pattern and its interannual variation, with the North Atlantic 
Oscillation as the first mode. EA/WR’ s climate impact extends from eastern North America to 
Eurasia. The positive (negative) EA/WR produces positive (negative) temperature anomalies 
over the eastern US, western Europe and Russia east of Caspian Sea, with negative (positive) 
anomalies over eastern Canada, eastern Europe including Ural Mountains and the Middle East. 
These anomalies are largely explained by lower- tropospheric temperature advections. Positive 
(negative) precipitation anomalies are found over the mid-latitude Atlantic and central Russia 
around ~60°E, where lower-level cyclonic (anticyc Ionic) circulation anomaly is dominant. The 
eastern Canada and the western Europe are characterized by negative (positive) precipitation 
anomalies. 

The EA/WR is found to be closely associated with Rossby wave propagation. Wave activity 
fluxes show that it is strongly tied to large-scale stationary waves. Furthermore, a stationary 
wave model (SWM) forced with vorticity transients in the mid-latitude Atlantic (~40°N) or 
diabatic heat source over the subtropical Atlantic near the Caribbean Sea produces well- 
organized EA/WR-like wave patterns, respectively. Sensitivity tests with the SWM indicate 
improvement in the simulation of the EA/WR when the mean state is modified to have a positive 
NAO component that enhances upper-level westerlies between 40-60°N. 
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1 . Introduction 


Large-scale teleconnection patterns that persist from days to months are known to play an 
important role in determining whether a particular season will be warm or cold, wet or dry. The 
teleconnection patterns, often defined in terms of the height and/or sea level pressure fields, and 
the corresponding atmospheric circulation fields typically cover vast geographical areas 
(Bamston and Livezey 1987). 

The North Atlantic Oscillation (NAO) is one of the primary modes of atmospheric variation 
over the North Atlantic, impacting seasonal climates over North America and Europe (e.g., 
Hurrell 1995; Donat et al. 2010; Sugimoto and Hanawa 2010). Other teleconnection patterns 
with an important component in the North Atlantic region are the Scandinavia (SC A) (Wallace 
and Gutzler 1981; Bueh and Nakamura 2007), the East- Atlantic (EA) (Washington et al. 2000; 
Bojariu and Reverdin 2002), and the East Atlantic / West Russia (EA/WR) (Bamston and 
Livezey 1987; Washington et al. 2000; Wang et al. 2011) patterns. While numerous earlier 
studies investigated the NAO because of its dominant large-scale influences on climate 
variability, impact of the other teleconnections on climate variability over the North Atlantic and 
generation of their patterns have yet to be addressed in more detail. Particularly, the EA/WR, 
which is characterized by two main large-scale anomalies located over the Caspian Sea and 
western Europe, has not attracted any significant and detailed investigation despite the 
recognition of EA/WR as one of the leading teleconnections over the North Atlantic (Barnston 
and Livezey 1987; Washington et al. 2000). Wang et al. (2011) and Lim and Kim (2013) 
suggested that understanding EA/WR is very important because not only its impact extends 
across the European mainland, but also the impact reaches mid-latitude East Asia as a planetary- 
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scale stationary wave pattern. The possible role of EA/WR in modulating the east Asian winter 
monsoon variability was analyzed by Wang et al. (2011) and Kim et al. (2013). 

Regarding the interrelationship among the teleconnections over the Atlantic, Scherrer et al. 
(2006) investigated the possible relationship between the NAO and large-scale atmospheric 
blocking located over the Atlantic, the Scandinavian Peninsula, and mainland Europe. 
Specifically, the study identified that a change in the intensity of the westerlies over the central 
latitudes of the eastern North Atlantic and over much of Europe associated with the NAO, is at 
times associated with a long-lived atmospheric blocking in the vicinity of Great Britain and the 
Scandinavian Peninsula. Pavan et al. (2000) and Shabbar et al. (2001) addressed that atmospheric 
blocking over the North Atlantic seems more frequently observed during a substantial negative 
NAO. However, there are some indications that blocking-like persistent high-pressure systems 
over the European mainland are more related to the positive NAO phase (Wanner et al. 2001). 
Overall, it appears that there is still much to be learned about possible inter-relationships 
between the various modes, including the potential role of other key large-scale patterns (e.g., 
ENSO) in modulating teleconnection patterns in this region. 

Turning to the physical mechanisms responsible for these patterns of variability, there is 
some evidence for an important role of wave propagation forced by transient eddies. Bueh and 
Nakamura (2007) found that transient eddies migrating along the storm track together with 
incoming Rossby wave activity from the Atlantic are crucial for maintaining teleconnection 
pattern. Several other studies also argued that heating over the Atlantic may be responsible for 
producing the atmospheric teleconnection (e.g., Walter et al. 2001; Wang et al. 2011). Clearly, 
the relative importance of heating and weather transients, and the primary source regions for 
generating Rossby waves associated with the various Atlantic teleconnections needs clarification. 
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The purpose of this study is to provide a better (physically-based) understanding of the 
impacts of the EA/WR teleconnection on winter climate over the North Atlantic and surrounding 
regions, and their relation to Rossby wave propagation. This study takes advantage of the latest 
high-resolution (-0.5° latitude/longitude) reanalysis data to document the impacts of the EA/WR 
pattern on temperature and precipitation over the Arctic, western Russia, eastern North America, 
North Africa and Europe. The stationary wave model of Ting and Yu (1998) is used to examine 
the role of stationary Rossby waves in explaining the structure and forcing of these patterns. 

Section 2 describes the data and stationary wave model utilized in this study. The capture of 
the EA/WR spatial pattern and its interannual variation over the study domain is described in 
Section 3. Section 3 also discusses the atmospheric circulation and other features associated 
with the EA/WR that acts to produce the regional climate impacts. Section 4 identifies the source 
regions and the propagation of large-scale stationary waves associated with the EA/WR pattern, 
and includes an assessment of the sensitivity of the response to the basic state tied to variations 
in the NAO or ENSO. This is followed by the concluding remarks in Section 5. 

2. Data and Model 

This study uses the Modem-Era Retrospective analysis for Research and Applications 
(MERRA) reanalysis data (Rienecker et al. 2011) to analyze the past 33 winters (December, 
January, February or DJF) from 1979/80 through 2011/12. The key variables consist of SST, 
upper-level (250hPa and 300hPa) geopotential height, wind (300hPa and 850hPa), temperature 
(300hPa, 850hPa, and 2 meter level), sea level pressure (SLP), diabatic heating (residually 
diagnosed), and precipitation. The horizontal resolution is 0.5° latitudex0.6667° longitude, and 
the temporal resolution is daily though the daily values are averaged to create monthly means. 
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The stationary wave model is a fully nonlinear baroclinic model with 14 vertical levels on 
sigma coordinates with R30 truncation in the horizontal (Ting and Yu 1998). The model 
variables include vorticity, divergence, vertical velocity (sigma coordinate), surface pressure, 
geopotential height, and temperature. A rigid-lid boundary condition is applied at the top and the 
surface of the model atmosphere. For damping, Rayleigh friction and Newtonian cooling are 
applied in the vorticity, divergence, and temperature equations to ensure meaningful solutions. 
For further details of the model see Ting and Yu (1998). 

3. Impact of the EA/WR pattern on temperature and precipitation 

The leading modes of large-scale teleconnections are computed from the monthly mean 
upper-level (250hPa) geopotential height field. The climatologies for each month of DJF were 
first removed from the raw data. A rotated EOF (REOF) technique (Richman 1986) was applied 
to the resulting anomaly data spanning the past 33-winters. The first four REOFs, in order of 
decreasing variance, are the NAO (Wallace and Gutzler 1981; Barnston and Livezey 1987), the 
SCA (Washington et al. 2000; Bueh and Nakamura 2007), the EA (Barnston and Livezey 1987; 
Bojariu and Reverdin 2002), and the EA/WR (Barnston and Livezey 1987; Washington et al. 
2000). Please note that a previous study using an alternative rotation method also captured the 
NAO as the first mode followed by SCA (Hannachi et al. 2009). The eigenvectors (the spatial 
distributions) are shown in the left panel of Figure 1, with the corresponding principal 
component (PC) time series shown in the right panel. Note that all spatial patterns and PC time 
series are plotted for what is conventionally considered to be the positive phase of the pattern. 
Also shown in the right panel are the teleconnection indices archived at the National Oceanic and 
Atmospheric Administration (NOAA)/National Center for Environmental Prediction 
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(NCEP)/Climate Prediction Center (CPC) (ftp://ftp.cpc.ncep.noaa.gov/wd52dg/ 
data/indices/t ele index.nh). The PC time series and the index time series are in good agreement, 
confirming that the previously identified teleconnection patterns have been successfully captured. 
This study has examined the sensitivity of the patterns to the domain size (e.g., to include the 
entire northern hemisphere) and found that the spatial distributions of the REOFs are robust with 
respect to changes in the domain (figure not shown). Note that an AO (Thompson and Wallace, 
1998) mode was captured as the fourth REOF (not shown), though in our regional analysis, the 
NAO likely already includes some of the variability of the AO. 

The EA/WR pattern (REOF 4), the main focus of the present study, has two large-scale 
anomaly centers located just north of the Caspian Sea and western Europe, influencing Eurasian 
climates (Figure Id). During the positive phase of the EA/WR, a negative height anomaly 
occurs over the Atlantic near 40°W and 40-45°N, while a positive anomaly is found over central 
Europe (0-30°E). The EA/WR appears to originate in the North Atlantic and extends 
northeastward across Europe and European Russia (Barnston and Livezey 1987; Washington et 
al. 2000). 

In order to quantify the temperature and precipitation anomalies associated with each 
teleconnection pattern, and to help explain how those anomalies are produced in different regions, 
several atmospheric variables at low levels (i.e., near surface and 850hPa) are regressed onto the 

time series associated with each z250 hPa height PC. In particular, the regressed field )) 

for the Mth mode at grid point (x,y) is defined as 

R M (.x,y) = Y.t= 1 T(x,y, t) ■ P M (t ), (1) 

where T (x, y, t)T(x,y ) is the anomaly field in question at time step t, and Pm (t) represents the 
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normalized monthly PC time series of the z250 hPa height for the Mth mode. In the above 
summation, nt is equal to 99 months, the length of the analysis period. 

In this study, the advective temperature change arising from the circulation anomaly is 
calculated as ~^Tel ■ Vr c; . (Linkin and Nigam 2008), where Vjei indicates the horizontal winds 
associated with the EA/WR pattern and Ten refers to the climatological temperatures (Figure 2c). 
Note that the contribution from the nonlinear component, “^.i ' VTrei , and the other linear 
component, -Vcu ' V7Y #l , is by comparison small. 

In general, the geographical distribution of the variables in Figure 2 reflects the upper-level 
height anomaly distribution that has a wave structure spanning the Atlantic, Europe, and western 
Russia shown in Figure 1 (cf., Barnston and Fivezey 1987; Washington et al. 2000). Such a wave 
structure is expected to be associated with circulation anomalies and temperature advection 
patterns that produce regions of alternating warm and cold weather. Figure 2b clearly shows a 
strong anticyclonic circulation over Europe and cyclonic circulation over western Russia in the 
event of positive EA/WR. These anomaly pattern helps to determine temperature advection 
characterized by warm advection over the northeast Atlantic extending into western Europe and 
western Africa (~20°W~10°E), and Russia northeast of the Caspian Sea. Cold advection occurs 
primarily over eastern Canada, and Russia near the Ural Mountains (Figure 2c), where northerly 
flow located west of the negative SLP anomaly is dominated. 

The resulting temperature anomaly pattern in Figure 2 a reflects this spatial distribution of 
circulation and temperature advection with the strongest response occurring over Russia. Figure 
2a reveals that the 2-meter air temperature (T2m) anomalies associated with the EA/WR extend 
to North America, the Middle East, Africa, and Eurasia. Western Russia shows larger magnitude 
of T2m anomalies (>1°C) than the magnitude of anomalies over the other regions (Figure 2a). 
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During the positive EA/WR, positive anomalies are found over the eastern US, western Europe, 
and Russia east of Caspian Sea. Eastern Canada, far eastern Europe including Ural Mountains 
(~60°E), the Middle East and northeastern Africa regions are characterized by negative 
anomalies, although the anomalies are not statistically significant at the 10% level over part of 
the regions. 

Figure 2d shows that distribution of precipitation anomalies is strongly coupled with the 
anomalous SLP and lower-level circulation pattern (Figure 2b). The positive precipitation 
anomalies are clearly found where negative SLP anomalies are located, and the opposite is true 
for the negative precipitation anomalies. During the positive EA/WR, the mid-latitude Atlantic 
and the central Russia (~60°N) are the main regions that have the above-average precipitation, 
while the eastern Canada and European region experience the below-average precipitation. A 
weak positive precipitation anomaly is observed over the eastern US, but the magnitude is 
relatively small. The anomalies described in Figure 2 are of course by definition of opposite sign 
for the negative phase of the EA/WR. 

4. The Connection with Rossby Waves 

a. Wave activity fluxes 

In this section the relationship between the EA/WR teleconnection pattern and large-scale 
Rossby waves is examined. Wave activity flux (WAF) vectors are calculated to identify the 
spatial distribution of wave propagation associated with the EA/WR. Following Plumb (1985), 
the stationary WAF is given as 
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where the variables ( u , v), p, T, and ® represent the zonal and meridional wind, pressure, 
temperature, and geopotential height, respectively, and where X and (p represent longitude and 
latitude. The constant co is the earth’s rotation rate (=7.292 x 10 5 rad s’ 1 ) and a. is the radius of 
the earth. The prime denotes the deviation from the zonal mean at each latitude and height. 
S = y- + y- is the static stability; the caret indicates an areal average over the Northern 

Hemisphere; k (=287 J K’ 1 kg‘Vl004 J K’ 1 kg’ 1 ) is the ratio of gas constant to specific heat at 
constant pressure, and H is a constant scale height. 

Figure 3 shows the horizontal WAF distribution at 300hPa. The pattern suggests that the 
EA/WR reflects Rossby waves responding to forcing in the mid-latitude Atlantic. Specifically, 
there is a clear wave train extending across the Atlantic, western Europe and Russia associated 
with the EA/WR. 


b. Rossby wave source 

The source regions and the propagation characteristics of the EA/WR teleconnection are 
investigated in this section using the stationary wave model (SWM) of Ting and Yu (1998) (see 
section 2 for the model description). This study focuses on both diabatic heat forcing and 
transient vorticity forcing, the importance of which in generating large-scale wave trains in the 
extratropics has been shown by previous studies (Sardeshmukh and Hoskin s 1988; Qin and 
Robinson 1993). Also, the upper-level WAF pattern (Figure 3) suggests that the Rossby wave 
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source (RWS) is located over the extratropical Atlantic, in a region where transient eddies 
associated with the Atlantic storm track could play a role (Bueh and Nakamura 2007). Previous 
studies emphasize the importance of changes in wintertime extratropical transient eddy activity 
for forcing changes in upper-level jet intensity and its extension (e.g., Losada et al. 2007). 

Following Schubert et al. (2011) this study attempts to reproduce the observed 
teleconnection patterns by forcing the SWM model with a series of regional heat and vorticity 
forcing functions, respectively, located at 5° longitude/latitude intervals distributed throughout 
the North Atlantic. The 3-dimensional basic state for this run is the MERRA climatology 
computed from the past 33 winters (1979/80-2011/12). The similarity is assessed between the 
SWM response to each forcing function and the observed teleconnection by calculating their 
spatial correlation over the North Atlantic and Europe domain (60°W-60°E ,20-80°N). These 
correlations are then plotted at the forcing locations (Figure 4). 

Figure 4a shows that the largest correlations for the SWM response to heat forcing are found 
over the subtropical western Atlantic near the Caribbean Sea. Correlation values over the mid- 
latitude Atlantic are generally small. The subtropical region where the largest correlation is 
found is consistent with the region of the largest positive diabatic heating anomaly (Figure 4b) at 
mid-troposphere and SST anomaly (Figure 4c) regressed onto the EA/WR, indicating a possible 
heat source over the region to generate the Rossby-type wave that resembles the EA/WR. This 
possible positive relationship between SST anomaly over the subtropical western Atlantic and 
the phase of EA/WR is in good agreement with Wang et al. (2011). The study suggested that 
cold North Atlantic SST anomalies are associated with the negative EA/WR pattern that 
enhances the Siberian high and east Asian winter monsoon. 
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As for the transient vorticity forcing, Figure 5 a shows that the EA/WR has spatial 
correlations greater than 0.7 over the mid-latitude Atlantic basin (~40°W, ~40°N), suggesting 
they have a substantial Rossby wave component that resembles the EA/WR pattern driven by 
extratropical transients. It also suggests that this wave pattern is most easily forced in the central 
North Atlantic near the climatological jet exit region (please see the superimposed contours in 
Figure 5a). Blackburn and Hoskins (2001) found that the exit region of the Atlantic jet stream 
was important for cyclone growth and dynamically forced ascent. Distribution of upper-level 
(250hPa) daily zonal wind variance in Figure 5b shows the largest variance over this mid-latitude 
North Atlantic, indicating the strongest transient eddy activity around the wave source region. 

Figure 6 shows some examples of the SWM response to idealized forcing chosen to 
resemble the observed EA/WR. The location of the heat forcing and vorticity forcing, 
respectively, was chosen based on the location of the largest correlation values in Figures 4 and 5 
which is 80-70°W, 20-25°N for the heat forcing, and 45-35°W, 35-40°N for the vorticity forcing 
(Figures 6a, b). The SWM responses (upper-panels) and the observed patterns captured by 
REOFs (lower-panels) show substantial similarities. In particular, positive heat forcing in the 
subtropical North Atlantic (20-25°N) and vorticity forcing in the mid-latitude Atlantic (35-40°N), 
respectively, generates a stationary wave with maximum amplitude over the high latitude North 
Atlantic, the Scandinavian peninsula, and the Eurasian continent (Figures 6a, b), similar to the 
observed EA/WR pattern in Figure 6c. The SWM, however, also generates a low latitude 
response (Figure 6a) that is not found in the observed pattern. Overall, Figure 6 supports the idea 
that EA/WR pattern is associated with large-scale Rossby wave responses to diabatic heat 
forcing and vorticity forcing over the extratropical Atlantic. 
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Next, rather than deducing the RWS based on the SWM responses to idealized forcing, the 
RWS is estimated directly based on observation data (MERRA data). Following Sardeshmukh 
and Hoskins (1988) and Qin and Robinson (1993), the RWS is derived from the quasi- 
geostrophic vorticity equation as 

RWS= -v x ■ V(C + /) - 0T + f ) v ■ v x , (1) 

where V y is the divergent (irrotational) wind vector, < the relative vorticity, and f the Coriolis 
parameter. The linearized fonn of RWS can be written as 

RWSl=-R; • V(? + 0 - (l + f)V • V x - V • V x - V x • V(', (2) 

where the overbar denotes the climatological mean and the prime denotes the anomaly associated 
with the teleconnection pattern (e.g., EA/WR). The first and fourth term on the right hand side 
of (2) are associated with vorticity advection, whereas the second and third terms involve the 
generation of wave vorticity by the divergence of the divergent wind (i.e., vorticity stretching). 
Sardeshmukh and Hoskins (1988) and Qin and Robinson (1993) found that the first and fourth 
terms capture the main features of the tropical source whereas the second and third terms 
detennine the main features of the extratropical RWS (Seo and Son 2012). Figure 7a shows the 
distribution of the RWSs characterized by positive peak over the mid-latitude Atlantic. Location 
of the peak values at 50-30°W and 35-50°N is consistent with the positive spatial correlation 
maximum shown in Figure 5 a. 

In order to confirm that the regions forcing the EA/WR wave patterns are indeed in the 
locations we discussed in Figures 4 through 6, the EA/WR pattern is reproduced with the RWS 
forcing in the region outlined in Figures 7a. The region is upstream of the EA/WR pattern, and 
located where the magnitude of spatial correlation between the patterns of observed 
teleconnection and idealized model result (by transient vorticity forcing) is high (Figure 5a). The 
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responses to the forcing in those regions (Figure 7b) show a reasonable similarity with the 
observed EA/WR pattern captured by the REOFs (Figure 6c). 

c. Sensitivity to background flow change (NAO and ENSO) 

Since variations in the NAO impact the upper-level westerlies over the North Atlantic, this 
study examines whether the NAO has an impact on the other patterns through NAO-related 
changes in the base state. Table 1 shows the variance of the PCs of the EA/WR REOFs for 
different NAO phases. The results show that the variances of the PCs are larger on average in 
the event of a positive NAO. The dependency of the intensity of the EA/WR pattern upon the 
NAO phase indicates that strong wind and stonniness over the European mainland is more likely 
to occur during the positive NAO, as suggested by Donat et al. (2010). 

In order to more clearly demonstrate that the EA/WR pattern is indeed sensitive to the NAO 
phase, the SWM response to vorticity sources is re-examined with modified base states. In 
particular, this examination considers base states that have added to it +/- one standard deviation 
of the NAO pattern. Additional experiments were conducted by incorporating +/- one standard 
deviation of the ENSO component 1 . Figure 8 shows the changes made to the base state, with one 
standard deviation of the positive (negative) NAO (Fig. 8a) acting to accelerate (decelerate) the 
westerlies between 40-60°N over the Atlantic, and decelerate (accelerate) them between 20-40°N. 
El Nino also plays a role in displacing the Atlantic storm track and modifying the westerlies, but 
with smaller magnitude compared with the NAO impact (Figure 8b). 

Four SWM experiments were conducted to quantify the impact of the NAO and ENSO on 

the EA/WR. Figure 9 shows the EA/WR response to the positive NAO (Figure 9a) and negative 

1 Note that the ENSO component is captured as a higher REOF (not shown) occurring just after 
the EA/WR mode. 
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NAO (Figure 9b), respectively. The transient vorticity forcing is identical to that referred to in 
Figure 5. An EA/WR-like pattern is generated for both basic states (Figures 9a, b), though there 
are differences in the intensities of the ridges/troughs between the two wave patterns. The wave 
path is also different. In general, EA/WR is better organized and strengthened over the Atlantic 
and western Europe when the mean state is modified with a positive NAO component. This is, to 
a great extent, consistent with Wanner et al. (2001) and Scherrer et al. (2006) who found a 
greater possibility of strong pressure patterns (teleconnection) and atmospheric blocking over the 
European mainland during the positive phase of the NAO. The wave track is also more realistic 
with the positive NAO added to the base state (Figure 9a). For the case of the negative NAO, 
there is a northward displacement of the anomalies making the response less like the observed 
EA/WR pattern (Figure 9b). Also, the pattern shows relatively smaller amplitude especially over 
the Atlantic and western Europe. It appears that the stronger upper-level westerly jet at 40-60°N 
over the Atlantic during the positive NAO may facilitate the eastward propagation of the EA/WR 
over the European continent, hi contrast, the main anomalies during the negative NAO tend to be 
situated further west over Greenland due to weakened upper-level westerlies between 40-60°N. 

The above experiments were repeated for the modified mean states reflecting the ENSO 
impact. The results show some differences in the EA/WR pattern for the El Nino and La Nina 
cases, but the differences are considerably less than those for the NAO shown in Figure 9, 
reflecting the smaller change in the upper-level jet compared to that associated with the NAO 
(figure not shown). This suggests that a very strong ENSO is necessary to noticeably modulate 
the EA/WR patterns, and that in general the NAO is a greater factor in modulating those patterns. 

5. Concluding remarks 
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This study has investigated the EA/WR, the important winter atmospheric teleconnection 
pattern generated over the North Atlantic region. The focus is on providing further insights into 
the spatial scope and physical mechanisms of the impacts, as well as the dynamical mechanisms 
responsible for the existence of EA/WR. 

It was found that the four leading REOFs of the upper tropospheric height field consist of, in 
order of decreasing variance, the NAO, the SCA, the EA, and the EA/WR. There have been few 
studies that have examined in any detail the impact of the EA/WR on regional temperatures and 
precipitation across Eurasia. Here the present study found that the positive (negative) EA/WR is 
associated with a strong wanning (cooling) over Russia east of the Caspian Sea whereas it is 
associated with cooling (warming) over the Ural Mountains of northern Russia. The strong 
temperature response over those regions is linked to a surface pressure anomaly near the Ural 
Mountains and associated strong atmospheric circulation and thermal advection anomalies. Other 
impacts linked to the positive (negative) phase of the EA/WR are a warming (cooling) over the 
eastern US and western Europe, and a cooling (wanning) over eastern Canada, far eastern 
Europe and the Middle East, though these tend to be weak. It was also found in this study that 
precipitation anomalies are strongly coupled with the anomalous SLP and lower-level circulation 
pattern. Positive (negative) precipitation anomalies are pronounced over the mid-latitude Atlantic 
and the central Russia (~60°N), while the negative (positive) anomalies are mainly distributed 
over the eastern Canada and European region during the positive (negative) EA/WR. The eastern 
US region exhibits the positive (negative) anomalies, but the magnitudes are generally small 
during the positive (negative) phase of the EA/WR. 

The second part of this study examined the relationship of the EA/WR teleconnection to 
large-scale Rossby waves. It was found, using a SWM, that the EA/WR pattern could be 


16 


354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 


reproduced by forcing the model with vorticity transients over the mid-latitude North Atlantic 
(35-40°N) and with diabatic heating over the subtropical western Atlantic (15-20°N). The results 
appear consistent with Bueh and Nakamura (2007) in that the transient eddies migrating along 
the Atlantic stonn track play an important role for emanation of the Rossby waves from the mid- 
latitude Atlantic. The RWS calculation based on observation provides supporting evidence that 
the RWS regions over the Atlantic are at ~40°W and ~40°N for the EA/WR. It was further 
shown that it is the vorticity transients near the Atlantic jet region (~40°N, -40° W) that are 
important for generating the EA/WR-like response. 

This study found that the EA/WR-like Rossby wave responses in the SWM show some 
sensitivity to the base state, with better-organized eastward propagation and enhanced 
ridge/trough anomalies over the North Atlantic and far western Europe occurring when the base 
state includes a positive NAO component. This appears to be due to the enhancement of the 
upper-level westerlies (40-60°N) associated with the positive NAO. Less organized EA/WR 
pattern is simulated when the base state includes a negative NAO component. This finding 
appears to support the argument by earlier studies (Pavan et al. 2000; Wanner et al. 2001: 
Scherrer et al. 2006) that addressed weaker pressure patterns and associated blocking over the 
European continent during the negative NAO. 

The modulation of the EA/WR pattern by other key large-scale modes of varaibility such as 
ENSO was also examined. The modulation of the teleconnection patterns by ENSO (El Nino 
versus La Nina) is less pronounced compared with the NAO impact, suggesting a greater role of 
NAO in modulating those teleconnection patterns. The result emphasizes the importance of 
upper-level westerly wind changes over the North Atlantic in modulating these teleconnection 
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patterns, suggesting that the variability in these regions consists of a complex interplay of the 
leading patterns. 

This study provides fundamental assessment of the EA/WR pattern in terms of 1) its climate 
impact on surface temperature and precipitation over vast area across North America, North 
Africa and Europe including western Russia, and 2) relation of the EA/WR to Rossby wave 
propagation. Investigation of the EA/WR should further be complemented by more detailed 
studies, including predictability of the generation of EA/WR pattern. Better identification of 
extra-tropical heating (or cooling) over the western subtropical Atlantic and upper-level transient 
eddy activity near the mid-latitude central Atlantic should be a key factor for reliable prediction 
of the EA/WR pattern. 
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Table 1 Variance of the principal components for EA/WR with respect to the phase of 
NAO. Numeric range in each cell represents the 90% confidence interval . 
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EA/WR 

NAO (+) (17 yrs) 

1.28 

0.57-2.00 

NAO 0(15 yrs) 

0.72 

0.30-1.14 
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REOFs of the DJF Z250 (1979DEC-2012FEB) 
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Figure 1 . The first four rotated empirical orthogonal functions (REOFs) of the monthly 250hPa 
height archived from MERRA reanalysis. The data set consists of data for 33 winters from 
1979/80 DJF through 2011/12 DJF. The left panel represents the distribution of non-nonnalized 
eigenvectors whereas the right panel the corresponding PC represents the time series (solid line). 
Dashed lines denote the teleconnection pattern indices time series archived at 
NOAA/NCEP/CPC. 
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Regressed T2m, UV850, SLR, T850 adv. and Prep onto EA/WR 


A) T2m C) T adv. (850) 




Figure 2. Distribution of a) 2 meter air temperature anomalies, b) sea level pressure and 850hPa 
circulation, c) the advective temperature change ( • VT cli ) by 850hPa atmospheric 

circulation, and d) precipitation associated with 1 standard deviation in the positive EA/WR PC 
based on a linear regression. For the calculation of temperature advection, Vjei denotes the 
horizontal winds regressed onto EA/WR and Ten represents the climatological temperatures. 
Shaded are the regions where the anomaly values are statistically significant at 10%. Wind 
vectors statistically significant at 10% level are plotted thick. 
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H300 8c WAF associated with EA/WR 
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Figure 3. Distribution of the wave activity fluxes and geopotential height anomalies at 300hPa 
associated with the positive phase of EA/WR. Geopotential height anomalies statistically 
significant at 10% level are shaded. 
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Figure 4. a) Spatial correlations between the observational teleconnection pattern captured by 
REOF and stationary wave propagation pattern produced by the stationary wave model. The 3- 
dimensional basic state in the SWM is the MERRA climatology computed from the past 33 
winters (1979/80-2011/12). In order to generate the stationary wave at each grid point, diabatic 
heat forcings with the maximum at mid-troposphere are given, respectively, at the grid points 
with 5 -degree longitude-latitude interval over the Atlantic. Geographical domain for spatial 
correlation calculation is 60°W~60°E and 20°N~80°N. Correlation values are plotted on the grid 
points where the diabatic heat forcing is given for generating the stationary wave in the model, b) 
represents the observed distribution of residually diagnosed diabatic heating anomalies at mid- 
troposphere (700-300hPa) regressed onto the EA/WR. c) represents the distribution of SST 
anomalies regressed onto the EA/WR. 
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B) U25Q variance 
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Figure 5. a) Same as Figure 4a except that transient eddy forcings of vorticities are given, 
respectively, at the grid points with 5 -degree longitude-latitude interval over the Atlantic. 
Correlation values are plotted on the grid points where the transient eddy forcing of vorticity is 
given for generating the stationary wave in the model. Contours in a) denote the climatological 
upper-level (300hPa) westerlies, b) represents the distribution of upper-level (250hPa) zonal 
wind variance regressed onto the EA/WR. 
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Stationery wave model response to heat & vorticity forcing (EA/WR) 

A) Heat forcing (8Q~7QW,2Q~25N) B) Vort. forcing (45~35W t 35~40N) 
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Figure 6. Upper panel: Streamfunction (divided by 10 6 ) of the simulated large-scale stationary 
wave propagation (a and b) similar to EA/WR pattern, a) represents the wave propagation forced 
by diabatic heat source, while b) by transient vorticity forcing. Geographical location of the 
forcing given in the model is 80~70°W, 25~25°N for the diabatic heat forcing, and 45~35°W, 
35~40°N for the transient vorticity forcing, respectively. The 3-dimensional basic state in the 
SWM is the MERRA climatology computed from the past 33 winters (1979/80-2011/12). Lower 
panel: Observational teleconnection patterns for c) EA/WR obtained from REOF. 


29 


622 

623 

624 

625 

626 

627 

628 

629 

630 

631 

632 

633 

634 

635 

636 

637 

638 

639 

640 

641 

642 

643 



Figure 7. Upper panel: Rossby wave source (RWS) distribution of the observational large-scale 
stationary wave associated with EA/WR pattern. Source region is represented by positive values 
(red shading). Lower panel: The EA/WR-like teleconnection patterns reproduced by stationary 
wave model with the given RWS in the boxed regions shown in upper-panel. Boxed region is 
where magnitudes of spatial correlations found in Figure 5a are predominantly high. 
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A) 1 atd. of positive NAO (U3QQ) 
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Figure 8. Upper panel (a): Distribution of one standard deviation of the positive NAO 
component in terms of upper-level (300hPa) westerly. Note that the climatology is the average of 
the upper-level westerlies over 1979/80-2011/12 DJF period. Lower panel (b): Same as a) but for 
the positive ENSO (i.e., El Nino) component 
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A) EA/WR, Istd. of NA0+ 
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Figure 9. Simulated EA/WR-like Rossby wave train embedded in the modified climatological 
basis state where one standard deviation of a) positive NAO is added and b) negative NAO is 
added. Wave propagation patterns are plotted in terms of streamfunction (divided by 10 6 ). 
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